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Abstract 

Three-step cascade decays into two invisible particles and two visible particles via two intermedi- 
ate on-shell particles develop cusped peak structures in several kinematic distributions. We study 
the basic properties of the cusps and endpoints in various distributions and demonstrate that the 
masses of the missing particles and the intermediate particles can be determined by the cusp and 
endpoint positions. Effects from realistic considerations such as finite decay widths, longitudinal 
boost of the parent particle, and spin correlations are shown to be under control for the processes 
illustrated. 
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I. INTRODUCTION 



At the energy frontier, the LHC experiments are taking us to an unprecedented terri- 
tory of the Tera-scale physics beyond the Standard Model (SM). At the cosmo frontier, we 
have entered an era of precision cosmology. With much progress made in the two frontiers, 
we have to admit that our understanding of the Universe is still far from being complete. 
According to the precise measurements of the cosmic microwave background fluctuations, 
such as WMAP about 95% component of the current universe has never been directly 
detected in the laboratory. The dominant component 72%) is dark energy that is re- 
sponsible for the accelerating expansion of the universe 2J. The second dominant (~ 23%) 
is cold dark matter (CDM), which is assumed to be in a form of nonrelativistic matter but 
cannot be explained within the SM. Albeit its extraordinary success in explaining current 
experimental data with incredibly high precision, the SM is regarded as an effective theory 
below a certain scale. For example, theoretical unnaturalness of the SM, dubbed as gauge 
hierarchy problem, suggests new physics beyond the SM at the TeV scale. Therefore, it is a 
very intriguing possibility that such CDM components may appear in new physics models. 

Indeed, some new particle physics models have answers for the astrophysical question 
about CDM. One of the most popular scenarios is a thermal production of weakly interacting 
massive particles (WIMP) [3]. In this scenario, a stable particle X had been once in thermal 
equilibrium in the early history of the universe, but got frozen out as its reaction rate 
became slower than the expansion of the universe. The stability of the CDM particle over 
cosmic time is often due to an unbroken parity symmetry or a discrete symmetry. Under 
such a symmetry, the SM particle fields are in the trivial representation while new particle 
fields are in some nontrivial representation. The decay of the lightest new particle into SM 
particles is prohibited. The current observation highly suggests that the CDM particle has 
its mass at the electroweak scale and its couplings with a size of weak interaction. Some 



popular models with WIMPs are supersymmetric models with R parity [4J, the universal 
extra dimension (UED) model with Kaluza-Klein (KK) parity [Si], and the littlest Higgs 
model with T parity Q]. 

This WIMP with an electroweak scale mass is likely to be produced at the LHC. In 
hadron colliders, such weakly interacting neutral particles can be identified only by missing 
transverse energy. The measurement of its mass is of crucial importance to reveal the 
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FIG. 1: The antler decay topology of a parity-even particle D into two missing particles {X\ and 
X 2 ) and two visible particles {a\ and a 2 ). 

identity of the CDM, but this is a very challenging task at the LHC because such invisible 
particles are pair-produced with combinatoric complications and with large errors especially 
in jet energy measurements. In the literature, many new ideas to measure the CDM mass 
have been proposed {7], such as endpoint methods jg|, polynomial methods [si, H], and M T2 



methods 



ll|-ll4| 



Recently, we have proposed a new approach to measure the missing particle mass by using 



the singular structures in the kinematic distributions of the antler decay [15l. Il6|. The antler 
decay is a resonant decay of a parity-even particle D into a pair of the missing particles 
(X\ and X 2 ) and a pair of SM visible particles (a\ and a 2 ) through two on-shell parity- 
odd intermediate particles (B\ and B 2 ), as depicted in Fig. [U We have studied two kinds of 
singular structures, an endpoint and a cusp. The positions of cusps and endpoints determine 
the masses of the missing particle as well as the intermediate particle, if the parent particle 
mass tub is known from other decay channels directly into two SM particles 1 . 

There are a few interesting merits of this method: (i) the positions of the cusp and 
endpoint are stable under the spin correlation effects since it is purely determined by the 
phase space; (ii) a cusp as a sharp and non-smooth 'peak is statistically more advantageous 
to search than endpoints, and more identifiable to observe than kinks; (iii) the simple con- 
figuration of outgoing particles can reduce combinatoric complication which is commonly 
troublesome in many missing particle mass measurement methods; (iv) the derived analytic 
functions for some kinematic distributions are very useful to reconstruct the mass parameters 



1 This is possible since the particle D has even parity. 



3 




FIG. 2: The cascade decay topology of a parity-even particle D into two missing particles {X\ 
and X2) and two visible particles (a\ and 02). 

by best-fitting. 



In this paper, as a companion of Ref. 16| . we focus on another decay topology with two 
visible particles and two missing particles from a parity-even particle D: cascade decays 
shown in Fig. [21 In this process, the parent particle D sequentially decays into two particles 
through three steps in series, finally ended up with a missing particle X2. There are two 
non-trivial types of three-step cascade decay, according to at which step the first missing 
particle X\ is produced. Type I and Type II cascade decays have different cusp and 
endpoint structures. Unlike the symmetric antler decay case with one kind of intermediate 
particle, the cascade decay involves two different intermediate particles. We thus need to 
fix one more unknown mass, which requires more independent observables. The study of 
the basic properties of cusp and endpoint in various kinematic distributions to determine 
the unknown masses for the three-step cascade decay is our main purpose. The cusp in 
the invariant mass distribution of Type I decay has been discussed in the context of new 



physics models with the CDM particle stabilized by Z3 symmetry 171 ] . 

The rest of the paper is organized as follows. In Sec. [IT], we categorize all possible kine- 
matic variables from the four-momenta of the two visible particles. Section II I II deals with 
the Type I cascade decay. We present the expressions of cusps and endpoints of various 
kinematic distributions in a common case where m ai = m a2 = and mx 1 = mx 2 - The 
functional form of the invariant mass distribution is also given. The general mass case is 
to be discussed in the Appendix. In Sec. IIVI we present the corresponding results for the 
Type II cascade decay. Section [V] is devoted to realistic considerations such as the finite 
widths of the intermediate particles, the longitudinal boost of the parent particle D, and 
the spin correlation. We then conclude in Sec. IVIl 
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II. KINEMATICS OF CASCADE DECAY TOPOLOGY WITH TWO MISSING 
PARTICLES 



We consider the four-body cascade decay of a heavy particle D through three con- 
secutive steps. The cascade decay resulting in a single missing particle and three visi- 
ble particles has been extensively studied in the literature. In the MSSM, a good ex- 
ample is the process of q — > qx 2 — > q£ n £ — > Q^fXi- m the UED model, we have 
Q(i) zWq L^£ n q -»■ B^£ f £ n q. Here £ n (£ f ) denotes the near (far) lepton with 
respect to the parent particle. In principle, three observable particles provide enough in- 
formation to determine all the unknown mass parameters involved [3, S]. However, there 
are some difficulties in extracting proper information, especially because of combinatoric 
complications. It is hard to distinguish £ n from £f. Furthermore, the parent particle D is 
to be pair-produced due to its odd parity (or nontrivial representation), and thus there is 
always another decay chain in the same event. 

Here we consider the three-step cascade decay with two missing particles. The parent 
particle D is of even parity and thus its single production is allowed. The final states 
are simply two visible particles (ai and a 2 ) with missing transverse energy. There is no 
combinatoric complication when forming the invariant mass of two visible particles. In 
addition, if the rest frame of D in the transverse direction can be determined, the individual 
transverse momenta of ai and a 2 in the frame also show kinematic singularities which have 
additional information on the mass parameters of the system. Note that this decent feature 
relies on the information of -D's transverse motion. 

The cascade decays of D — )■ a 1 a 2 X 1 X 2 can be classified according to at which step the 
first missing particle, say X±, is produced. Note that we fix the other missing particle {X 2 ) 
to be produced at the last step. If X\ is also from the last step, the final intermediate 
particle B is just missing and this decay is indistinguishable from a two step decay. We do 
not consider this case. Then, there are two non-trivial three-step cascade decays, as depicted 
in Fig. |2j In the Type I decay, X\ is from the second step. The parent particle D decays 
into a visible particle ai and a new particle C, followed by the decay of C into a missing 
particle Xi and a new particle B. Finally B decays into a visible particle a 2 and a missing 
particle X 2 . In the Type II decay, X\ is from the first step: D decays into CX±, followed 
by C ai B, and finally B — > a 2 X 2 . In the view point of two observable particles a\ and 
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a 2 , this Type II decay is a two-step cascade decay of a new heavy particle C. As shall be 
shown, there is no cusp structure in Lorentz-invariant distributions. 

It is useful to describe the kinematics of the three-step cascade decay in terms of the 
rapidity of individual massive particles or a combination of multiple particles: 

v (k) = ^- (1) 

'It r k y \ L ) 

ml 

where E { and rrii are the energy and mass of the particle (system) i in the rest frame of a 
particle (system) k. To avoid confusion, we adopt the following rapidity notations for the 
Type I and Type II decays: 





Type I Cascade 


Type II Cascade 


rapidity notation 


6 


d 



For the sake of simplicity, when the rapidity is defined in the rest frame of its mother particle, 
we omit the superscript specifying the reference frame. 

With the four-momenta k\ and k 2 of the two observable particles a\ and a 2 in the lab 
frame, respectively, we consider the following observables in three categories: 

• Lorentz invariant observables: the invariant mass of a\ and a 2 , 

m = y/{h + k 2 y . (2) 

• Longitudinal-boost invariant observables: 

— the magnitude of the transverse momentum of a visible particle i, 

Pn = |kf | , (3) 

— the magnitude of the transverse momentum of the a\-a 2 system, 

p T =|kf + k^|, (4) 

— the transverse mass of the a\-a 2 system, 

m T = ^Pr + m 2 . (5) 

• Non-invariant observables: 
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cosine of 0j, the angle of the visible particle a« in the cm. frame of a\ and a 2 , 
with respect to their cm. moving direction, 

kj aa) ■ k( g ) 
|kf a) ||k(°)| 



cos e = . (6) 



Here the bold-faced letter denotes the three- vector momentum, k = k\ + k 2 , and the super- 
script (D) and (aa) denote the D-rest frame and the cm. frame of ai and a 2 , respectively. 

As shall be shown, p?i and mx distributions show cusp structures if the mother particle 
D is produced at rest in the transverse direction. At a hadron collider, this is possible if D 
is singly produced. These additional cusp structure are very valuable to determine all the 
unknown masses. However strong QCD interactions always yield, for example, sizable initial 
state radiation, which causes transverse kick to the mother particle D: the cusps in pTi and 
rriT distributions can be smeared. Caution is required when drawing the consequences for 
the mass measurement from these cusps. In addition, the cos B distribution is defined in the 
rest-frame of the mother particle D. At a hadron collider, this is not observable. In what 
follows, we assume that the mother particle is produced at rest in the transverse direction. 

In general, the involved seven particles (D, C, B, a%, a 2 , X\, X 2 ) may have different 
masses. In many new physics models, the cascade decay processes of interest have massless 
visible particles and the same kind of invisible particles. For most of the presentation in the 
main text, therefore, we consider only the following case: 

m ai = m a2 = 0, m Xl = m X2 . (7) 

The result for the most general case with seven different masses is presented in the Appendix. 

III. TYPE I CASCADE DECAY 

As illustrated in Fig. 121 the Type I cascade decay is the decay of a parity-even particle 
D into two missing particles X\ and X2 and two visible particles a\ and a 2 through 

► C + ai (h), (8) 
C^B + Xx, 

B — > a 2 (k 2 ) + X 2 . 
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Here the particles D, C, at and a 2 are parity-even while the particles B, X\, and X 2 are 
parity-odd. In order to accommodate the Type I cascade decay, we need at least two heavy 
parity-even particles. 

One good example for this decay channel is in the universal extra dimension (UED) 



model 



181 ] . It is based on a single flat extra dimension of size R, compactified on an S\jZ 2 



orbifold. All the SM fields propagate freely in the whole five-dimensional spacetime, and 
each field has an infinite number of KK excited states. Since the KK parity is conserved, 
the lightest KK particle (LKP) with odd KK parity is stable and becomes a good candidate 



18 



19] 



for the CDM. Usually the first KK mode of the U(l)y gauge boson B^ is the LKP 
All the second KK states of the SM particles have even KK-parity and mass of ~ 2/R. 
Lower limit of 1/R > 300 GeV is set based on the combination of the constraints from the 
p parameter [20j , the electroweak precision tests 2l| , the muon g — 2 measurement 22] , the 
flavor changing neutral currents 23j, and direct search by DO group at the Tevatron 24]. 
The second KK modes are within the reach of the LHC. Possible Type I cascade decays 
are 



g {2) -» q n + q {2) ^q n + B {1) -> q n + B {1) + q f B {1) . 



(9) 
(10) 



Now we present the cusps and endpoints of m, tut, Pt, Pn, and cos 9 distributions in 
terms of the masses. As in Eq. 0, two missing particles are of the same kind and the visible 
particles are massless in this case. The rapidities of the particles B and C in the rest frame 
of their mother particles are given by 



cosh £ 



B 



m c 
2m 1 



1 + 



m 



B 



rn 



x 



lb \ TTVq TJIq J Allbc \ UL D, 

We will also use E n and Ef, the energy of the near a\ and the far a 2 in its mother's rest 
frame, respectively: 



cosh £, 



c 



m D 
2uif 



1 + 



m, 



c 



m, 



(11) 



E„ 



m D 



1 _ a^c 



D 



E 



f 



m B 
2 



1 - 



m 



x 



m. 



(12) 



For illustration, we take three sets for the mass parameters in Table [B The Mass-Ai is 
motivated by the decay in Eq. ([9]). The KK masses are determined by the UED model 
parameters of AR = 20 and 1/R = 500 GeV, where A is the cutoff scale ljj . The equal 
spacing of the KK mode spectrum in flat extra dimension leads to very degenerate masses, 



S 





m D 


m c 


m B 


m x 


6? 


Mass- 


-Ai 


1045.7 


1023. 


514.2 


500.9 


0.12 


Mass- 


Bi 


600 


400 


200 


100 


0.60 


Mass- 


"Ci 


600 


500 


150 


100 


1.16 



TABLE I: Test mass spectrum sets for the Type I cascade decay. All masses are in units of GeV. 

i.e., rriD mc ~ 2m# ~ 2mj. The Mass-Bi has substantial gaps for each pair of adjacent 
masses. Finally Mass-Ci has a sizable mass gap between mc and rriB + mx- 

For precise mass measurements using the singularities, it is required to have a visible 
cusp and/or endpoint in a kinematic distribution. The visibility of the endpoints can be 
determined by the functional behavior near the endpoint. In what follows, the shape of an 
endpoint is to be classified into a fast dropping one and a long-tailed one. 
(i) Invariant mass m distribution: We first discuss the distribution of the invariant mass m 
of two visible particles. The differential decay rate dT/dm is 



dT_ 

dm 



^ / max\2 



for < m < , 
for m casT <m< m%£i, 



(13) 



where the cusp and endpoint are 



(m 



(14) 



Note that the functional behavior of dT/dm is the same as that of the antler decay 16]. 
More general case with 7 different masses is discussed in the Appendix. 

Whether this cusp is sharp enough to probe can be easily deduced from Eq. f JT3|) . The 
dT/dm function is linear in m for m < rn^f, and a concave function for m™gf < m < m^f. 
At m = m^i/e, the concave function reaches its maximum. If m^^/e < m™^, which is 
equivalent to < 1, the cusp can be considered to be pronounced. 

In Fig. |3j we show the normalized differential decay rate of dT/dm. In order to compare 
the cusp shapes only, we present it as a function of m/m™^. The vertical lines denote the 
positions of in units of m 1 ^^. The Mass-Ax case with = 0.12 has a very sharp 

m cusp. The Mass-Bi case with = 0.60 shows a triangular shape with a cusped peak. 
However, the Mass— Ci case with = 1.16 has a dull cusp. In the dT/dm distribution, the 
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0.2 0.4 0.6 0.8 1 

FIG. 3: The normalized differential decay rate of the invariant mass of two visible particles, 
for the Type I cascade decay. The masses are in Tabled! 
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FIG. 4: The normalized differential decay rate of the transverse mass of two visible particles, 
Pg^- for the Type I cascade decay. The masses are in Tabled 

profile shape near the endpoint can be generally regarded as fast-dropping, as suggested by 
Eq. (US]). 

(ii) Transverse mass mx distribution: Figure H] shows the rate of the transverse mass rriT 
distributions. For all three Mass— Ai, Mass-Bi and Mass-Ci cases, the distributions 
show visible cusp structures. It is interesting to note that the Mass-Ci case has a more 
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FIG. 5: The normalized differential decay rate of the transverse momentum of two visible particles, 
v d J pT for the Type I cascade decay for the masses in Table [IJ 

visible cusp compared with the m cusp. We also note that this is contrasted to the antler 
decay case where there is no cusp in the ttit distribution 16] . As shall be shown in the next 
section, the Type II cascade decay also has a cusp in the distribution. Therefore, the 
presence of the tut cusp can be used for identifying the cascade decay topology The cusp 
and maximum positions in terms of the masses are 



(m T )Zl = E n + E f <*>~<* 



casl 



(15) 



where E n and Ej are in Eq. (fl2|) . 

(Hi) The system pt distribution: Figure [5] shows the normalized distribution dT/dpx of the 
transverse momentum pt of two visible particles. For all three mass spectra in Table [H 
the pt distribution has smooth peak without a cusp structure. Still the endpoint of pt 
distribution can be observed, which is 



(p 



max 
T /casl 



E n + E f e~^ cHB . 



(16) 



Only the Mass-Ax case has a fast dropping endpoint shape, which is attributed to very small 
momentum transfer to the visible particles. More general cases of Mass-B! and Mass- 
Ci have long-tailed endpoints. The pt distribution is not useful for the mass measurement. 
(iv) Single particle pxi distribution: The individual transverse momentum pTi distributions 
of two visible particles show unique functional behaviors, as shown in Fig. [6j The thin solid 
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FIG. 6: The normalized differential decay rate of the transverse momentum of one visible particle, 
r ^p T . for the Type I cascade decay. In the Mass-Ai case, the line labeled by "near" ("far") denotes 
the pTi distribution of a\ (02). Thick lines are the summed distributions of pr%- 

line labeled by "near" ( "far" ) is the pTi distribution of the near visible particle a\ (the far 
visible particle a<i). The pxf distribution has both the cusp and the endpoint structures, 
while the pt u distribution has only an endpoint. This pT n endpoint has a sudden ending 
shape like a step function, which holds true for all mass cases. 

In most realistic situations, one may not distinguish the near visible particle from the 
far one. Here we show a more practical observable, the transverse momentum of any visible 
particle, which becomes the sum of both pxi distributions. The thick lines in Fig. [6]represent 
the sum 

The position of the cusp and the endpoint in the pxf distribution is given by 



and the endpoint in the pt u distribution is located at 



(17) 



/ \ max 
{PTn) casl 



(18) 



where E n , Ef, £ B and are in Eq. ffTl| and f fT2|) . Depending on whether (pjy,)™ > 
(Pt/)™si (^ ne case °f Mass-Ax) or not (the cases of MASS-Bx and MASS-Ci), the summed 
distribution shows apparently different shape, as shown in Fig. El While (pTn can be 
easily determined due to the unique spiky feature of the distribution, the cusp and the 
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FIG. 7: The summed distributions of cos 0j in the Type I cascade decay for the masses in Table 
I 

endpoint of the pxf distribution when (pTn)^a < (^ T /)casi * s ra ther difficult to identify in 
the summed pxi distribution. 

(v) cos0 distribution: The variable cosG in Eq. ((6]) is defined by the angle of one visible 
particle. We have two cos0 distributions for a\ and 02, which cannot be distinguished. In 
Fig. [7J therefore, we present the summation of two cos 0, distributions in the rest frame of 
D for Mass-Ai and Mass-Bi cases. It is symmetric about cosG = 0, and has two cusp 
structures, cos 6^1 1 an d cos O^asi 2 > marked by the vertical arrows. In terms of masses, they 
are 

E n - E f exp(-£ B - £c 



cos 



n cuspi = E n - E f exp(£ B -£c) 
casl E n + E f exp(Z B -&) 



E n + E f exp(-£ B - fa) 



In the Mass-Ai case, cos 0^1 1 stands on a steep slope, which is difficult to probe. The 
Mass-B! case shows two pronounced cusps. 



IV. TYPE II CASCADE DECAY 



Type II cascade decay is a chain decay of 

D{P)- 



C + X x , 
C — > B + ai (h), 

B — > a 2 (k 2 ) + X 2 . 



(20) 



13 
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m c 


m B 


m x 




Mass-a 2 


614 


299 


222 


161 


138.0 


Mass-b 2 


600 


300 


200 


100 


193.6 


Mass-c 2 


400 


250 


150 


120 


120.0 



TABLE II: Test mass spectrum sets for the Type II cascade decay. All masses are in units of 
GeV. 
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FIG. 8: The normalized differential decay rate of the invariant mass of two visible particles, 
for the Type II cascade decay. The mass spectrum sets are described in Table ILT1 

A good example can be found in the MSSM: 

H/A^xi+Xi xi^L + i, i^tf + Xy (21) 

As in the Type I cascade decay, we restrict ourselves to the realistic cascade decay with 
m ai = m a2 = and mx 1 = mx 2 - Then there are two independent rapidities, (c and (b- 

coshC S = ^ l + -r « coshCc = ^ U + -T--T • 22 
2m b V m c J 2m c V m D m D J 

For illustration, we take three mass sets for the Type II cascade decay in Table [III 

(%) Invariant mass m distribution: We first study the distribution of the invariant mass of 

a\ and a 2 . Note that in the view point of a\ and a 2 , this Type II cascade decay is a three 

body decay of the parent particle C. The presence of the invisible Xi decayed from D does 

not change any Lorentz invariant result. The m distribution is the same as that of, i.e., mm 
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of the decay X2 ~^ ^ — > ^fXi m the MSSM. This m« distribution is well known to have 



no cusp structure. The endpoint is 



m. 



cas2 / 
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m. 
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x 



(23) 

b cJ \ m bJ 

In Fig. [BJ we show the m distribution for three sets of the mass parameters in Table [Til a U 
of which have right-angled triangle shoe without a cusp. 

The absence of a cusp in a two-step cascade decay can be understood by a simple kine- 
matic configuration. For the antler decay (D — > B\ + B2 — > a>iX\ + 02X2) in the massless 
visible particle case (m ai = m„„ = 0), the following four critical points correspond to a 



kinematic singular structure 161 ] : 
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(24) 



Here we simplify the picture as an one-dimensional case. It is clear to see that happens 
when two observable particles move in the same direction, while two kinematic configurations 
of back-to-back moving correspond to either or m™ s ^ 2 . For a two-step cascade decay 

(C — > a± + B — > ai + 02X2), ai and a<i in one- dimensional space have only two independent 
kinematic configurations, moving in the same direction and moving in the opposite direction. 
There is no critical point left for the cusp. 

(ii) Transverse mass ttit distribution: Unlike the invariant mass distribution, the distri- 
bution contains the information about the transverse momenta of the first missing particle 
X\. As shown in Fig.[9j there is a cusp here. We stress once again that this cusp appears 
only when D is produced at rest in the transverse direction. 

Another interesting feature is that the position of the cusp is nothing but the m 
maximum: 



fm \ CUS P — m max 



(25) 



This non-trivial equality is a unique feature of the Type II cascade decay. 
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FIG. 9: The normalized differential decay rate of the transverse mass of two visible particles, d mT 
for the Type II cascade decay. The mass spectrum sets are described in Table HH 
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FIG. 10: The normalized differential decay rate of the transverse momentum of two visible 
particles, r dp T f° r the Type II cascade decay. The mass spectrum sets are in Table ITTl 

(in) System distribution: The total pt distributions for the Type II cascade decay are 
shown in Fig. [TUJ All three mass sets have smooth pt distributions. And their endpoints 
are all long-tailed. This feature is common for the antler, Type I, and Type II cascade 
decay topology. 

(iv) Single particle pTi distribution: Figure HT1 shows the distribution of the individual trans- 
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FIG. 11: The normalized differential decay rate of the transverse mass of one visible particle, 
Yd^ T . for the Type II cascade decay. The left figure is for the near visible particle a%, and the 
right one is for the far visible particle a2- 



verse momentum of the near a\ and the far a%. The near p^ n distribution has a sharp cusp 
and a fast dropping endpoint. However the pxj distribution has a long tailed endpoint 
without any cusp. In terms of masses they are simply 



Note that the product of (pm)cas2 an d (.PTn)™2 removes the Cc* dependence, which depends 
on the intermediate masses mc and mg. In addition the ratio (prn)^! / (Prn)SS depends 
only on the rapidity (c- 

As discussed before, the individual pTi distribution cannot be constructed. Instead we 
show the sum of two distributions in Fig. [T2J For the Mass-a 2 case, the cusp in the pT n 
distribution and the smooth peak of the p?f distribution are located nearby. In their sum, 
the pm cusp survives over the relatively round prf peak and the fast dropping pT n endpoint 
is also measurable. For the Mass— C2 case, however, the pm cusp and the pr/ peak are 
separated so that the summed distribution shows both. With finite number of data, it 
would be difficult to distinguish the px n cusp from the pxj peak. 

(v) cos0 distribution: We consider the cosG distribution for the Type II cascade decay. In 
Fig. [T3]we show the normalized dT / 'dcosQ for the near and far visible particles (denoted by 
thin lines) as well as their sum (thick lines) for the Mass-A 2 and MASS-C2. In both cases, 
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FIG. 13: The sum of dT/d cos 0j for the Type II cascade decay. 



the summed distribution of cos Oj is symmetric about cos = 0, and has one independent 
sharp cusp denoted by vertical lines in Fig. [13J The cos cusp position in terms of the mass 
parameters is 



cos cusp 

Lub w cas2 



m c 1 - 



m c 1 



+ m B 1 



(27) 



(vi) Mass determination from the cusps and endpoints: Unlike the antler decay with one 
kind of intermediate particles, the cascade decay has two different intermediate particles. 
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In addition, the Type II decay has fewer independent observables of cusps and endpoints. 
First there is no m cusp structure. Second the rriT cusp position is the same as the m 
endpoint. A natural concern arises whether we have enough information to determine all 
the masses, especially at the LHC where the cos cusp cannot be used. Fortunately three 
unknown masses {mc, mg, and mx) are unambiguously determined by three observables of 
™£5, (Pt„)™ S 2 P , and (pmK 



where R n is 



and «i 2 3 are 



m c = R a m D , m B = x jl- — m C) rn x = Jl- — m B , (28) 



«3 — «1 — «2 



I max\2 

= ^cas 2 ) (30) 

1 ^> /7 ; / \ cusp ' v 1 

'cas2 



Q 2 



2m gv /(p rn ) c m a l /(p Tn ); 

V^Tn)^/^ 



cusp 
TnJcas2 



o > , (PTn)S | / (PTn)cas2 

(PTn )cas2 V (PTn)^S 

V. EFFECTS OF REALISTIC CONSIDERATIONS 

All the previous expressions of the cusps and endpoints have been derived in an idealistic 
situation: the total decay widths of decaying particles are ignored; the D rest frame is 
assumed to be reconstructed; the spin-correlation effects from the full matrix elements are 
negligible. In this section, we investigate these effects on the position and shape of each 
kinematic cusp and endpoint. 

A. Finite width effects 

Up to now we have applied the narrow width approximation, ignoring the width of de- 



1511 . we focus on the effects of 



caying particles. Since the effect of finite is very minor 
T B and T c . 

We find that the mass spectrum is the most crucial factor to determine the stability of 
the cusp and endpoint structures under the width effects. Out of six cases in Tables HI and HTl 
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FIG. 14: The finite decay width effects on the m and mx distributions in the Mass— Ai case. 
Solid lines are for the on-shell decay, the long dashed lines for T/M = 0.01, the short dashed lines 
for T/M = 0.1, and the dotted lines for T/M = 0.5. Here T/M = T B /M B = T c /M c . 




FIG. 15: The finite decay width effects on the summed distributions of pTi and cos0, in the 
Mass-Ai case of the Type I cascade decay. As before, we take T/M = 0, 0.01, 0.1, 0.5. 

the Mass-Ax has very vulnerable structures. This case is special because of its degenerate 
masses: the observable particles have very small momentum transfer and their kinematic 
phase space is highly limited. 

In Fig. [TH and Fig. [151 we show the finite width effects for the Mass-Ax case. We present 
four cases for T B and T C - on-shell (solid line), T/M = 0.01 (long dashed line), T/M = 0.1 
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(short dashed lien), and T/M = 0.5 (dotted line). Here T/M = To/mB = Tc/rric for 
simplicity. Just one percent of T/M destroys all the sharp cusp structures into smooth 
peaks. In addition, the positions of the peaks are shifted significantly from the true cusp 
positions. There is no way to extract the mass information from the cusps. For T/M > 0.1 
the summed pTi and cos Gj distributions lose their functional behaviors completely, leaving 
very smooth and featureless distributions. 

The fast-falling endpoints in the m, rriT, and pxi distributions are also smeared out con- 
siderably. The degree of its shifting is large even for T/M = 1%. One interesting observation 
is that two shifted endpoints of the m and distributions are the same to be - 2mx, 
denoted by vertical arrows. This new endpoint is from the kinematic configuration where 
two visible particles' momenta span all the phase space determined by mo and rrix- Even 
though we do not know the intermediate particle masses, the missing particle mass mx can 
be read off. For this information, the distribution is more advantageous than the m 
distribution, because of its fast falling shape. 

In a realistic new physics process, however, this Mass-Ai case does not allow even one 
percent of T/M. For example, the decay in the mUED model has the decay widths of 

T D = T z(2) ~ 270 MeV, T c = T m ~ 5 MeV, T B = T lW ~ 1 MeV, (31) 

which leads to T/M ~ 10~ 5 . This is attributed to the limited phase space. In summary, 
the extreme Mass— Ax case has generically negligible width effects. The cusp and endpoint 
structures are reserved. 

We consider more general mass parameters, Mass-Bx for the Type I and MASS-C2 for 
the Type II cascade decay. First we examine the finite width effects on the invariant 
mass distributions in Fig. [Tj3 These cases show more stable cusp and endpoint structures 
from the finite width effects. For T/M = 1%, the m distributions in both Type I and 
Type II decays do not change, keeping the same cusp and endpoint structures. For 10% 
of T/M, the m cusp of the Type I decay retains its position, though losing its sharpness. 
The m endpoints in both Type I and Type II decays are shifted into the new position 
mo — 2mx- If T/M = 50%, the Type I decay does not retain the shape and position of 
the cusp, and the Type II decay does not show the right-angled triangle shape of the m 
distribution. Both cases have the same new endpoint at mo - 2mx, which is also valuable 
information for the missing particle mass measurement. 
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FIG. 16: Finite width effects on the normalized m distribution. We take the Mass— Bi case for 
Type I decay, and Mass-C 2 for Type II decay. As before, we take T/M = 0, 0.01, 0.1, 0.5. 
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FIG. 17: The width effects on the normalized dT/dmx for the Type I Mass-Bi and 
Type II Mass-C 2 cascade decays. As before, we take T/M = 0, 0.01, 0.1, 0.5. 



In Fig. [T71 we show the width effects on the rrij- distributions. The % cusp structures 
are more stable than the m cusps in both Type I and Type II decays. For T/M = 1%, 
the changes in the distribution are unnoticeable. For T/M > 10%, we start to lose the 
sharpness of the cusps but still keep the positions for the cusp in both cases. If T/M = 50%, 
the cusped peaks become dull further with relatively stable positions, and the my endpoints 
are shifted into mp — 2mx- 
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FIG. 18: The width effects on the summed distributions of pxi for T/M = 0, 0.01, 0.1, 0.5. 
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FIG. 19: The summed cos 6 distributions for r/M = 0, 0.01, 0.1, 0.5. 

Figures [18] and [19] show the width effects on the summed distributions of p?i and cos 6j 
respectively. Both distributions preserve the cusp structure for T/M = 1%. If T/M > 10%, 
however, the finite width effects almost smear the cusp and endpoint structures. 



B. Longitudinal boost effect 

In hadronic collisions, the longitudinal motion of the particle D is not determined, which 
affects only the cos0. The angle Oj of the visible particle aj is defined in the cm. frame of 
a\ and with respect to their cm. moving direction, and this direction is defined in the D 
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FIG. 20: Normalized differential decay rates versus cos© in the L>-rest frame (thin curves) and 
in the pp lab frame with y/s = 14 TeV (thick curves). 

rest frame. 

In order to see the longitudinal boost effects, we convert the cos O distribution in the 
D rest frame into the pp frame at the LHC, by convoluting with the parton distribution 
functions of a proton. In Fig. [201 we compare the summed distributions of cos 0j in the 
.D-rest frame (thin curves) with that in the pp lab frame at \fs = 14 TeV (thick curves). 
For the parton distribution function, we have used CTEQ6 26|. We take the Mass-Ax for 
Type I and the MASS-C2 for Type II decay. For simplicity we assume that the heavy 
particle D is singly produced through the s-channel gluon fusion and qq annihilation. 

Unlike the finite width effects, the longitudinal boost effect does not completely smash the 
characteristic shape. The sharp cusp structures survive to some extent in both Type I and 
Type II cascade decays. The shift of the cos cusp position is minor. Moreover the overall 
functional shape remains the same even though the absence of events around cos = in 
the D rest frame is filled by the longitudinal boost effects. The cusp in the cos distribution, 
though Lorentz non-invariant, is quite useful to draw mass information. Again we emphasize 
that the e + e~ linear collider does not have this unambiguity. 
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with and without spin correlations. 

C. Spin-correlation effect 



Our main results are based on the kinematics only, ignoring the spin-correlation in the full 
matrix elements. Since this paper is focused on the basic properties of the kinematic singular 
structures in the cascade decays, full analysis for each new physics process is beyond the 
scope of this paper. Nevertheless the algebraic singularity origin of the cusp and endpoint 
keeps them stable under the spin correlation effects j^. 

In order to demonstrate this, we consider one example, the decay in the the UED 
model: 



Cascade Type I: Z {2) ^£ + L {2) ^£ + -> ££B^B m . 



(32) 



In Fig. EU we show their spin correlation effects. We found that the spin correlations do not 
change the m and pTi distributions. Two distributions with and without spin-correlation 
effects are almost identical. 
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Antler 


Cascade 
Type I Type II 


m 


yes 


yes no 


rriT 


no 


yes yes 


PT 


no 


no no 


PTn 
PTf 


yes 


no yes 
yes no 


cos G 


yes 


yes yes 



TABLE III: The presence or absence of the cusp in the kinematic distributions of m, mr, Pt, PTi, 
cos of the antler, Type I cascade, Type II cascade decays. 

VI. SUMMARY AND CONCLUSIONS 

We have studied the singularity structure, such as cusps and endpoints, in the kinematic 
distributions of three-step cascade decay of a new parity-even particle D and the determi- 
nation of the missing particle mass by using such singularities. 

Two non-trivial decay topologies, called the Type I and Type II cascade decays, have 
been studied. In the Type I decay (D — > aiC, C — > X±B, B — > a 2 X 2 ), where the first 
missing particle X\ is from the second-step decay, the distribution of the invariant mass m of 
two visible particles, a\ and a 2 , develops a cusp. Full functional form of the m distribution 
for general mass parameters has been derived. If the mother particle D is produced at rest 
in the transverse direction, various longitudinal-boost invariant observables accommodate 
cusp structures. First there is a cusp in the transverse mass rriT distribution, which is 
complementary for the m cusp since the rriT cusp shape is sharp even when the m cusp 
is dull. Although the transverse momentum distribution of the cm. system of two visible 
particles does not develop a visible cusp structure and a sharp endpoint, we note that the 
transverse momentum distribution of the far visible particle a 2 has a cusp, and that of the 
near visible particle a\ has an endpoint of the shape of a steep cliff. We also study the 
summed distribution of cos6j, which has two independent cusp structures. 

In the Type II decay (D — > X X C, C — > aiB, B — > a 2 X 2 ), the first missing particle Xi is 
from the first step decay. The kinematics of the two visible a x and a 2 is determined solely by 
the two-step cascade decay from the first intermediate particle C, so that the invariant mass 
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distribution does not have a cusp structure. However, the kinematic distributions regarding 
the transverse motion from production of both X\ and X2 can carry the information from 
the whole three-step cascade decay. We show both the % and YliPTi distributions have 
distinctive cusp structures. In the individual transverse momentum distribution, only the 
near visible particle has both a sharp cusp and a fast-falling endpoint. The cos G distribution 
also shows a cusp as well. Including the antler decay topology, we have summarized the 
existence of cusp in the kinematic distributions of m, tut, Pt-, Pn, cosG in Table [mi 

We have also considered the effects of finite decay widths, longitudinal-boost of the parent 
particle D, and spin correlation. The effects of the finite widths of the intermediate particles 
can be significant if the decay width is sizable, say Y/M > 10%. As the decay width 
increases, the sharp cusp gets smeared, and the endpoint position gets shifted to m D — 2mx- 
the missing particle mass mx can be still extracted by a proper fitting. The longitudinal 
motion of the parent particle D affects the distribution of cosG. At least for the sample 
mass parameters, however, the cos G cusp remains sharp after convoluting with the parton 
distribution functions of a proton at the LHC. Spin correlation effects from full S'-matrix 
elements turn out to be negligible in most cases, which is expected since the singularities 
are determined by the kinematic relations. 

With the companion paper on the detailed study of the kinematic cusps in the antler 
decay [if]], our analysis shows the general kinematic properties and provides useful formulae 
for the decay topologies with two visible particles and two missing particles. By looking 
at the singularity structures of various kinematic distributions, the hidden nature of the 
missing particle can be probed effectively and elegantly. With the outstanding performance 
of the LHC and detectors, this is an exciting time for such investigation. 
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Appendix: Invariant mass distributions for the general Type I case 



In this appendix, we present the invariant mass distribution in the general Type I cascade 
decays: 



»• C + oi(A;i), 
C— ►fl + Xt, 

B — > a 2 (k 2 ) 



(A.l) 



As discussed in the main text, the Type II cascade decay is practically a three body decay 
in the view point of visible particles. This four-body decay has generally seven different 
mass parameters. We define the rapidities of six particles as 



cosh £ c 
cosh £b 
cosh£ a2 



m 2 D 


+ m 2 c — 






2m£,mc 


■ 


m 2 c 


+ m 2 B - 


m 2 Xi 




2mcmB 




m% 


+ m l 2 - 


m 2 x 2 



2m B m 



cosh£ ai = 
cosh£ Xl : 
cosh £x 2 



m 2 D + m 2 



m 



c 



2m D m ai 


■ 


m 2 c + m 2 Xi — 


m% 


2m c m Xl 


? 


m 2 B + m 2 X2 - 


m l 2 



(A.2) 



a 2 



2m B m X2 



A very useful kinematic variable is \i the rapidity of the particle a 2 in the rest frame of a\\ 

11 i 
i , m — mi — mi 

= t^z (A-3) 



X 



2m ai m a2 



where the superscript (a\) denotes that the rapidity is defined in the rest frame of a\. 

The functional expression of dT /dm is different according to the mass relations. The 
derivation of dT /dm is similar to that presented in the appendix of Ref. 16|. For simple 
presentation, we introduce 



i++ = + £ai + ia 2 + Cc, 

£+- = \£b + £ai — £a 2 — £c\, 

£-+ = \Cb — Cai + Ca 2 + £c|> 

£, — ICfl — Coi — ia 2 — £c|- 



(A.4) 
(A.5) 
(A.6) 
(A.7) 



We order £_ + and £ and name them £i < £ 2 < £3- Analytic functions forms of dT /d\ 

are then written as 
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dT 

dx 



oc < 



if |6? - 62 - 6*1 > 61 Or £ B + 62 + 6* < 61 , 

— 6 + cosh -1 x, if cosh^i < x < cosh£ 2 , 

6-6, if cosh £ 2 < X < cosh 6, 

£ ++ - cosh -1 x, if cosh ^ 3 <x< cosh£ ++ , 

0, otherwise. 

if 16? - 62 - fci < 61 < 6; + 62 + fc, 

2cosh _1 x, ifl<X<cosh6, 

6 + cosh -1 x, if cosh 6 < x < cosh£ 2 , 

6 + 6, if cosh 6 < X < cosh 6, 

£++ - cosh -1 x, if cosh 6 < X < cosh 

0, otherwise. 



dT 

dx 



oc < 



(A.8) 



(A.9) 



The positions of the minimum, cusp, and maximum of the invariant mass distribution 



are 



M mm 
Ly± casl 



m ai + m l 2 + 1rn ai m a2 cosh 6, for H lt ... >6 
m ai +m a2 , for 7^7,... ,12 



m ai + m a 2 + 2m ai m a2 cosh ^3, 



y/ + + 2m ai m a2 cosh£ ++ . 



(A.10) 
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